Hall effect in the perovskite manganites 
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We compute the zero temperature phase diagram and the Hall response of the doped per- 
ovskite manganites within the model of double exchange and Jahn- Teller coupling, employ- 
ing the d — ► oo approximation proposed by Millis et ai. We find that in this scenario a 
"hole - like" Rh for the hole doped manganites, as observed at low temperature by Matl et 
alB, can be obtained only when lattice distortions persist in the metallic state as T — > 0. 
The calculated temperature dependence of Rh seems to be consistent with the measured 
"normal" Hall coefficient in thin films. 
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I. INTRODUCTION 



With the current interest in "colossal magnetoresis- 
tance" the ferromagnetic manganese oxides have been 
intensely studied over the last few yearsu. These com- 
pounds, of the form La\- x {Sr,Ca) x MnO^, are antifer- 
romagnetic insulators in the undoped phase [x = 0), and 
have a ferromagnetic ground state beyond a critical dop- 
ing. The ferromagnetism is conventionally understood 
in terms of the double exchange modeou. The transi- 
tion from the ferromagnetic to the paramagnetic phase 
is accompanied by a simultaneous metal-insulator 'tran- 
sition'. The transport is typically activated in the para- 
magnetic phase. 

The large increase in resistivity near T c and the rather, 
low transition temperature prompted the suggestion^ 
that, in addition to the double exchange mechanism, 
carrier localisation due to electron-phonon coupling ef- 
fects could -be relevant in these materials. Detailed 
calculationsffij incorporating both double exchange and 
the Jahn- Teller coupling have been successful in repro- 
ducing most of the qualitative features of the data. So, 
although there are possibly complications due to electron 
correlation effects as well as disorder due to cation sub- 
stitution, the basic interplay between charge, spin and 
lattice degrees of freedom in the manganites seems to 
be adequately described by the double exchange/ J aim- 
Teller scenario. 

Experimentally there is overwhelming evidenced for 
dynamical lattice distortions in these compounds, includ- 
ing isotope effects on T c and direct measurement of the 
Debye- Waller factor. There is also a systematic and large 
change in lattice constant with temperature. The acti- 
vated transport in the paramagnetic phase has been in- 
terpreted in detail in terms of small polaron hopping. 
While it is obvious and acknowledged that the electron 
phonon coupling is central to understanding the "high 
temperature" , T > T c , physics, its effect on the low tem- 
perature metallic phase is much less studied. This spin 
polarised phase is believed to admit an almost "band 
theoretic" description, with electron-magnon interactions 
describing the transports. The measurement of the Hall 
effectB in Lao.67Cao.33-/VfnO3 at low temperature, how- 



ever, suggests that the electron-phonon coupling may be 
crucial in understanding this phase too, as we discuss 
next. 



II. HALL MEASUREMENTS ON THE 
MANGANITES. 



The Hall response of the manganites has signatures 
of various interaction effects. There is strong field de- 
pendence to the Hall resistivity, as is usual in ferro- 
magnets, and the data can be parametrised as p xy = 
R H (T)B + R S (T)M, where R H and R s are respectively 
the "normal" and "anomalous" Hall coefficients, B the 
induction and M the spontaneous or induced magneti- 
sation. Measurements indicate that Rh > 0, while the 
anomalous term Rs < Q_ /in the Ca doped systemu) and 
\Rs\ > Rh for T ~ TJB- In fact, at all but the lowest 
temperatures, p xy is dominated by the anomalous term. 
The sign and temperature dependence of Rs is intrigu- 
ing and it has been suggested that its origin may lie in 
a combination of spin=£>rbit coupling effects and strong 
electronic correlationstil. Here we concentrate on the nor- 
mal Hall coefficient which, although much smaller than 
i?s, also has several mysterious features, not necessar- 
ily linked to the ferromagnetism; (a). Rh is "hole- like" 
although a simple electron count with degenerate Mn or- 
bitals would predict an "electron- like" Hall coefficient. 
(6). The carrier density, inferred from R H , at x = 0.33 
seems to be five times the nominal doping. So, either the 
Drude form RJj 1 ~ x, naively expected for the doped in- 
sulator, is inapplicable, or there is a serious discrepancy 
between the nominal and actual carrier density. There- 
fore, beyond the anomalies attributable to spin orbit cou- 
pling effects, the doping and temperature dependence of 
the normal Hall coefficient itself in the manganites is not 
understood. We have calculated Rh within the model of 
double-exchange and Jahn- Teller coupling introduced by 
Millis et al. as a function of electron-phonon coupling, g, 
electron density, n, and temperature, T. 
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III. MODEL 

We use the standard Hamiltonan, 

H = H e i + Hd-ex + H p h + H e i-ph 

The electronic band is spec- 

ified by H e i = J2ij t"f 4a«7 C 30<r ~ Y,iV n i' where a and 
(3 refer to the two e g orbitals. The double-exchange term 
Hd-ex = —Jh Yli S l c -<Ji couples the conduction electrons 
to the core 'spin' 3/2, which we treat classically, assum- 
ing JnSc/t — > oo. The classical phonons are described 
by a harmonic distortion term H p h = X), /2, and the 
Jahn- Teller coupling H e i^ ph = g £V Ti.c ioi<T Ta.pCip a where 
t is xt x + yr y , and r's are Pauli matrices. This model, 
even with spins and phonons treated classically, can be 
solved exactly only in tho_eJ — > oo limit. We follow the 
prescription of Millis et alS3 in handling the problem and 
the reader should consult their papers for the detailed 
formulation within the d — > oo framework. We skip the 
technical details here, referring the interested reader to 
our appendix for a discussion of the T = limit of the 
equations. 

For following our results it is sufficient to know that 
this scheme provides us with the distribution of lat- 
tice distortions and local spin orientation, P(r,6), and 
the electronic Greens functions G arT (k,uj), for electrons 
with spin parallel/anti parallel to the magnetisation. 
The transport coefficients are computed from the spec- 
tral functions, A aa (h,u>) = — (l/Tr)ImG (T(T (k ) u>) asllia, 
o xx = &H + °~ll an( l °xy = o'H + &l y where, for ex- 
ample, all = / depo{e) J duj(—d f / 'du>)A^(e^,Lj) , and 
°U = I depa{e)e J du(-df/doj)A^(e^u), p Q {e) is the 
free band density of states, and / is the Fermi function. 
The low field Hall coefficient, as usual, is Rh — &xy/oxx- 
We work in units where the half-bandwidth, D — 2. 

Unlike in [12] where most of the calculations were done 
at particle-hole symmetry, we have to work away from 
(rii) = 1 to obtain a non trivial Hall coefficient. We 
have solved the self-consistent equations numerically for 
various g, n and T and discuss the results next. 



IV. RESULTS 

A. T = 0, phase diagram, and doping dependence of 

Rh 

At low temperature, T -C T c , the spin background is 
almost completely polarised. In the absence of any or- 
bital ordering, the phonon distribution depends only on 
the scalar distortion, r = \f\, and, for T — > 0, the prob- 
ability distribution, P(r), for the (classical) phonons is 
sharply peaked, since the variance (Ar 2 ) - T. There- 
fore, in this limit, the electronic spectrum is completely 
determined by the average onsite distortion (r) = vq. 



This permits a simple classification of the various regimes 
as T — > 0. Following [12] these are (a). Fermi liquid; 
(r) — > as T — > 0, electron-phonon interaction effects 
vanish and there are no "quenched" distortions, (b). 
"Rigid band" phase, < (r) < r c (n), where r c (n) de- 
notes the critical distortion beyond which the system, 
with electron density n, becomes a polaronic insulator, 
and (c). (r) > r c (n) corresponding to the insulating 
state. Obviously, since phononic effects completely dis- 
appear in the first case, the Hall effect would simply cor- 
respond to that of a spin polarised band, i.e the band 
theoretic result. In the third case, since the conductivity 
vanishes as T — > 0, the Hall coefficient is ill-defined. The 
most interesting regime corresponds to (6), where lattice 
distortions persist in the T — > metallic state and Rh 
can significantly differ, in sign and magnitude, from the 
band-theoretic value. The metallic or insulating charac- 
ter of the ground state is determined from the density of 
states, N(uS), at co = 0. For the Fermi liquid this cor- 
responds to the band DOS, for the 'rigid band' phase it 
is finite but suppressed from the band value, and for the 
polaronic insulator N(Q) = 0. 

Our result for the "phase boundaries" separating the 
Fermi liquid, the distorted metallic phase and the pola- 
ronic insulator at T — are shown in Fig.l. Before dis- 
cussing the behaviour of Rh in the various phases let us 
establish the band theoretic result. Consider half-filling 
first, ignoring the antiferromagnetism for the moment. If 
the two e g orbitals were equally occupied, and the spin 
states were degenerate, then the "half-filled" state would 
correspond to quarter filling of spin degenerate bands, 
and the Hall effect would be "electron- like" . However, if 
the T — > state were spin polarised the majority carriers 
would inhabit two orbitally degenerate half-filled bands. 
This state, for a particle-hole symmetric bandstructure 
would correspond to Rh = 0. Hole doping on this state, 
i.e a reduction of the electron count, would lead to an 
electron-like Hall effect, not hole-like as has been mea- 
sured. 

The numerically obtained T — > result for Rh con- 
firms that for g < g c (n = 1) ~ 1.14, where the half- 
filled state becomes insulating, the Hall coefficient, at all 
doping, correspond to the band prediction. These states 
all have vanishing distortion as T — > 0. The results on 
Rh, for several interaction strengths, g = 0.9, 1.0, 1.3, 
are superposed in Fig. 2. Between g c (l) and g' c (l) the 
half-filled state remains insulating while states away from 
half- filling are metallic, but with finite (r). However, for 
a given g, with g c (l) < g < g' c (1), beyond a critical 
deviation away from half-filling the distortion seems to 
disappear. In this region, bounded by g c (n) separating 
the Fermi liquid and distorted phase, Rh exhibits inter- 
esting deviation from the band value. Notice that while 
Rh = would occur only at the particle-hole symmetric 
point, n — 1, in the Fermi liquid phase, it occurs away 
from half-filling in the distorted phase. This is the key 
to the anomalous behaviour in the doped insulator. 

We have shown a typical plot, for g = 1.3, in Fig. 2, 
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to illustrate this behaviour. Notice that with approach 
to half-filling, where the interaction effects are strongest 
due to the lowest kinetic energy, Rh moves off the band 
result, changes sign somewhere, and tends to diverge as 
n — ► 1. It is numerically difficult to push the imaginary 
frequency calculation below T ~ 0.01, but an asymp- 
totic analytic solution (see Appendix) at T — indicates 
that for x — > 0, Rh ~ x^ 1 / 3 . Increasing g beyond g' c (i) 
states away from half-filling rapidly turn insulating. The 
boundary g' c (n) separating the polaronic insulator from 
the Fermi liquid and distorted metallic phases is shown 
in Fig.l. 



B. Temperature dependence of Rh 

Electron-phonon coupling leads to a suppression of 
electron kinetic energy, as does magnetic disorder within 
the double exchange model. So, qualitatively, the effect 
of increasing temperature, reducing magnetisation, and 
suppressing electron motion, would look like increasing g 
at T — 0. The finite T results bear this out, as shown 
in Fig. 3. For the chosen coupling constant, g = 1.5, 
for x — 0.2 the hole-like Rh at T = becomes larger 
as T — > T c , for x = 0.3 the low temperature negative 
Hall coefficient changes sign to hole-like at some inter- 
mediate temperature, while for x = 0.4 it becomes less 
negative with increasing T but never changes sign. A 
complementary situation would hold for 'electron-doped' 
systems, if they could be experimentally realised. Over- 
all, Rh deviates even more from its band theoretic value 
with increasing temperature (T < T c ). 

Within our approximation the magnetic disorder is 
maximal at T = T c (the spin distribution becomes 
isotropic) and there is no further suppression of kinetic 
energy for T > T c . In this regime the behaviour of Rjj 
is probably determined by the activation effects in a xx . 



V. DISCUSSION 

A 'hole-doped' insulator is intuitively expected to have 
a hole-like Rh, so our approach might seem to be an elab- 
orate way of arriving at an obvious conclusion. The state 
we describe is, however, rather non trivial, and below we 
contrast it both to (a) the doping of a simple band in- 
sulator and (b) the mean field picture of the "distorted" 
metallic phase (where every site has the same distortion 
r). Next, (c), we contrast the properties of this mean- 
field state to that with randomly oriented distortions and, 
finally, (d), comment on quantum effects which we have 
neglected in our calculation. 

(a). Usually "electron" or "hole" doping refers to the 
doped carriers occupying a specific part of the Brillouin 
zone. For a band insulator, the doped holes, by defi- 
nition, occupy the hole like parts of the zone, and lead 
to a positive Rh- When the insulating state arises out 



of interaction effects, at an electron density where band 
theory would have predicted a metal, it is not clear at all 
that a reduction in electron density is equivalent to "hole 
doping" in the earlier sense. Also, if Rh as T — » 0, 
away from the insulating state, is required to deviate 
from the simple band value, some signature of interac- 
tion effects must persist in the metallic ground state. We 
have stated this earlier in the context of our results, here 
we would like to emphasize that, quite independent of 
specific models, a change in the sign of the Rh at T = 
requires "quenched" distortions of some sort. Proximity 
to an insulator is necessary, but not sufficient, to obtain 
the behaviour observed. 

(b) . Since the Jahn- Teller coupling can lead to a lat- 
tice distortion and splitting of the electronic band, naive 
band theory predictions based on degenerate e g levels are 
incorrect, as we have argued before. One could, however, 
construct a mean- field picture incorporating Jahn-Teller 
distortions, which might well provide an .-adequate de- 
scription of the half-filled insulating statellil. By conti- 
nuity one might argue that a similar mean-field picture 
for Jahn-Teller distortions away from half-filling, with a 
self-consistently computed r would provide a reasonable 
description of the metallic phase. In that case one would 
indeed dope into the 'lower band', and the quasiparticles, 
in this orbitally ordered phase, will show a hole- like Rh- 
That situation, though conceptually appealing, does not 
seem to be experimentally relevant. There is no signa- 
ture for orbital long range order in the metallic state of 
the manganites and one is required to deal with spatially 
random distortions, as we have done. 

(c) . The phase we consider has some distinct, non 
mean-field like, properties. Within the scheme we have 
adopted for the transport calculation, the single parti- 
cle spectral function completely determines the trans- 
port coefficients. Since the spectral function A(k,u>) = 
— (\/ir)Im{uj + /i — e£ — S(w)} -1 , the low frequency be- 
haviour lim^^oifJ- — ~ ft ~ miT, determines the 
low temperature transport. Here p = fi — £_r(0) and 
r = — (l/7r)E/(0) These parameters, p, and T have a 
physical interpretation. T is the single particle scatter- 
ing rate, and determines the resistivity within this model, 
while e£ = p determines the location of the Fermi sur- 
face and the Hall coefficient. This distinction between the 
intuitive mean-field picture and the orbitally disordered 
phase shows up in the behaviour of these parameters, 
and have direct physical consequences. The implications 
are (i) When the electrons propagate in a background 
of random distortions, they will experience a large scat- 
tering rate, in contrast to the orbitally ordered phase 
where V and p(T) -» as T -» 0. (it) The orbitally 
ordered metallic phase will show a Fermi surface in pho- 
toemission, while the state computed by us does not have 
any discontinuity in its momentum distribution function 
since, near half-filling, since one finds p > 2D the band- 
width (see Appendix) . It is somewhat difficult to cleanly 
verify these effects in the electronic spectrum because 
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there is already large substitutional disorder, leading to 
a finite r as T — > 0, and a smearing of the Fermi surface. 

(d). We have made a classical approximation for the 
phonons and, even in the absence of orbital ordering, be- 
low some coherence scale the quantum effects could be- 
come important and anneal out the random lattice dis- 
tortions. In this regard we can gain some insight from 
what is known about doped Mott insulators. In the Hub- 
bard model at large U, over an intermediate temperature 
range, one can view the system as electrons interacting 
with static random local moments. This doped Mott in- 
sulator has a doping dependence of Rh similar to what 
we have found here. However, below T ~ xD, where x is 
the doping, a Fermi liquid ground state is recovered and 
Rh becomes band like again. A similar situation could 
arise in this model also, as we outline below. Suppose 
the relevant optical phonon in our case has a frequency 
lo = y/k/m. For T <C fiui^/kg we could integrate out 
the phonons and generate a term in the Hamiltonian of 
the form J2i( c ial c ip) 2 ■ At T = the model reduces 
to a "Hubbard" model involving two orbital species cL Q 

with an effective onsite repulsion U e f f = U cou i + (at 
finite T we need to consider the other spin species too 
since the ferromagnet is not fully polarised). This sug- 
gests that there could be a temperature scale below which 
our results, involving static distortions, would not hold. 
Apart from the presence of substitutional disorder, it is 
not clear why such a quantum scale is not visible in the 
manganites. 

We can now discuss how the measured Rh fits within 
this framework. Since experimentally the half-filled state 
is insulating the relevant regime for us would be g > 
<7 C (1). Let us choose g = 1.3, as shown in Fig. 2, for il- 
lustration. There Rh is indeed large and positive near 
half-filling. However, the behaviour changes for x ~ 0.3 
where Rh ~ and changes sign at larger doping. This 
is effectively the interpolation between a doped insulator 
(strong coupling, low kinetic energy) and 'band metal' 
(weak coupling, large kinetic energy). So the enhance- 
ment seercl in R^ 1 at x = 0.33 could be simply due to 
the proximity of Rh to a zero crossing. For our choice of 
g, x = 0.33 itself is electron like, but that depends on the 
exact value of the coupling, which is anyway unknown. 
The two predictions one can make within this scheme 
are; (a). With increasing doping, within a given chemi- 
cal family, Lai- x Ca x MnOz say, Rh should change sign 
from hole-like to electron like, unless some other order- 
ing phenomena intervenes (b) . With increasing coupling, 
Pb substitution on Get, keeping the doping level fixed, 
Rh should deviate further from its band-theoretic value 
(should become larger if it is already hole- like, etc).. 

On the T dependence of Rh itself, there does not 
seem to be any published data. Unpublished results from 
Bryan Lin et a/Jl3 for LaQ.^CaQ.^MnO^, indicates that 
the hole-like Rh at T — increases with T and has a peak 
near T c , beyond which it seems to fall off. This seems to 
be consistent with our calculated T dependence (x = 0.2, 



Fig. 3). However, pn measurements on the manganite are 
rather difficult, due to the large magnetoresistance, and 
it would be useful to have independent confirmation of 
this trend. 

We would like to thank Boris Shraiman, Y.B.Kim, Art 
Ramirez and Harold Hwang for several discussions. 

VI. APPENDIX; THE T -> EQUATIONS 

As we have seen, much of the physics of Rh{x,T) is 
controlled by the coupling constant and doping depen- 
dence at T = 0, so it may be helpful to consider the 
T = equations in some detail. We essentially follow 
the derivation in [12]. Since the spins are completely po- 
larised, it is sufficient to consider only the Greens func- 
tion for electron spin parallel to the magnetisation, Gf-f. 

Within the d — > oo framework one solves a local prob- 
lem, of a site hybridising with a self-consistently com- 
puted bath, with the local Green's function G 1oc (oj) being 
the local projection, Ga(ui), of the lattice Green's func- 
tion. Broadly, for our electron-phonon problem, there 
are three steps; (a). Compute the local Green's function, 
G|° c (r, oj) for a given value of local distortion, r. This is 
elementary since the local fermion problem is quadratic. 
(b). Compute the full local Greens function by averag- 
ing over coordinates, with the appropriate distribution 
function. 

G\°{{<jj) = J drP{r)G\°f{r,uj) 

This is possible since the phonons are classical, (c). Cal- 
culate the distribution function, P(r), as the electronic 
free energy for distortion rasP~ e -TrinG(r) _ ^phis c l ose s 
the self consistent loop. 

More specifically, for the semicircular band DOS that 
we have assumed, the following equations need to be 
solved for T -> 0. 

If 1 1 

a(iu) n ) = iLu n + n - - / drP(r){ 1 ■ } (1) 

2 J a- gr a + gr 

which combines steps (a) and (b) above. a(ui) is related 
to G\f(ui) by 

G${w) = \f drP(r){^— + — L_ } 
1 2 J a — gr a + gr 

The phonon distribution 

P(r) = -^-exp[~ + J2 H al ~ g2 f )] (2) 

Zloc 2T ^ {UjJn) 

and Zioc = J drP{r). A more complicated version of 
these equations hold at finite T, where the distribution of 
local spin orientations also has to be taken into account. 

To understand the doping dependence of transport co- 
efficients we need to solve the equations above in the case 
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where P(r), in the limit T — > 0, reduces to a delta func- 
tion at one non zero value, i.e, 



a(iu) n ) 



/' 



1 



gr 



gr ' 



-} 



(3) 



and r is determined by maximising P(r). We can view 
this as a cubic equation for a(cu), for a specified gr. It 
is possible to solve this equation near the band edge, 
parametrised in terms of the deviation of the chemical 
potential \x from \x* (where one begins to dope into the 
lower band), and the DOS looks like, 



-(iMlmG 1 ^) 



- (m* - m)0(w - 0* - m)) 



From this the number of doped holes can be estimated 
to be x = 1 - n ~ (/i* - ^) 3 ^ 2 - 

To obtain <7 xa; etc, we need to use the relation 

H - E(w) = G ret (w) + G^tH 

which, for i?eG(0) ^ 1, leads to JmS(O) ~ JmG(O) ~ 
(/i* — /i) 1 ^ 2 . Combining this with the result for n be- 
fore, we obtain T ~ x 1/3 . Also /j - i?eE(0) ~ ReG(0) + 
iteG _1 (0) > 2 since 7mG(0) -> 0. 

Since T — > while ^ — i?eE(0) — ^ 0, the spectral 
functions, which enter the transport calculation can be 
roughly approximated as A(k,ui) ~ r/(eg — /2) 2 . From 
which, employing the formulae cited previously in the 



text, one obtains o x 



r 2 



i.e p 



~ 2 / 3 , and ct T 



r 3 . 



This leads to Rh ~ cr xy /a 
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Figure Captions 

Fig.l: Phase diagram at T = with varying electron 
density, n, and Jahn- Teller coupling, g. The phase above 
the solid line, g' c (n), is a polaronic insulator. The shaded 
phase, bounded by g c (n) and g' c {n), has (r) ^ and is 
metallic except at half-filling where it is insulating. The 
dashed line corresponds to Rh = which occurs away 
from half-filling in the distorted metal. In the region 
bounded by Rh = and g' c {n) Rh has the 'wrong' sign 
with respect to band theory. 

Fig. 2: Density dependence of Rh (arb. units) at T = 
for various J-T coupling, g = 0.9 (circles), g = 1.0 
(squares), g = 1.3 (triangles). The coincident g = 0.9 
and 1.0 traces are the same as the band theory result 

Fig. 3: Temperature dependence of Rh (arb. units) for 
g = 1.5 and different dopings, x = 1 — n. Results for 
x = 0.2 (circles), 0.3 (squares) and 0.4 (triangles). 
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